Anisotropic Thermal Conduction in Supernova Remnants: Relevance 
to Hot Gas Filling Factors in the Magnetized ISM 
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ABSTRACT 

We explore the importance of anisotropic thermal conduction in the evolution of supernova remnants 
via numerical simulations. The mean temperature of the bubble of hot gas is decreased by a factor of 
^ 3 compared to simulations without thermal conduction, together with an increase in the mean density 
of hot gas by a similar factor. Thus, thermal conduction greatly reduces the volume of hot gas produced 
over the life of the remnant. This underscores the importance of thermal conduction in estimating the 
hot gas filling fraction and emissivities in high-stage ions in Galactic and proto-galactic ISMs. 

Subject headings: conduction — MHD — supernova remnants — X-rays: ISM — ISM: magnetic fields 



1. Introduction 

Supernovae play an important role in driving the ther- 
mal evolution of the interstellar medium (ISM hence- 
forth) . The interiors of supernova remnants (SNR hence- 
forth) are filled with extremely hot plasma, with temper- 
atures of 10^ — 10^ K. However, the interaction of these 
remnants with the ISM is not well understood. Obser- 
vational clues to the nature of this interaction are made 
possible by the spectacular maps and spectroscopy of x- 
ray emission provided by the Chandra space telescope 
(Hwang 2004) , and by X-ray spectra from Chandra and 
FUSE of high-stage ions like O vi (Bowen et al. 2005; 
Oegerle et al. 2005; Savage & Lehner 2006). 

The time evolution of the volume of the hot gas pro- 
duced by a SNR has important implications for the ther- 
mal evolution of the multiphase ISM, as described by 
Field et al. (1969), Cox & Smith (1974), and McKee & 
Ostriker (1977) (M077 henceforth). A comprehensive 
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Table 1: Initial conditions of the ISM in our simulations. 
p is the density in amu cm~'^; T is the temperature in 
K; B is the magnetic field in /xG; and L is the grid size 
in parsecs. 
Run p T B L 



LO 


0.7 


8000 


0.0 


300 


LI 


0.7 


8000 


3.0 


300 


HO 


5.0 


10000 


0.0 


200 


HI 


5.0 


10000 


6.0 


200 



study has been performed by Slavin & Cox (1992, 1993) 
(SC92, SC93 henceforth) in one dimension. M077 had 
estimated that the hot gas would occupy 60-70% of the 
volume of the ISM. The models of SC92 and SC93 found 
that the porosity of hot gas is reduced by thermal con- 
duction relative to the prediction of the M077 model, 
suggesting that the warm phase, and not the hot phase, 
might be the pervasive phase of the ISM. The volume 
filling factor as well as the amount of gas that shines 
in high-stage ions such as O vi and O VIII has recently 
been measured by FUSE and Chandra (Oegerle et al. 
2005; Yao & Wang 2005; Savage & Lehner 2006). Such 
data has made it possible to subject the M077 model to 
strong observational tests. Consistent with these tests, 
numerical modelling of the phase structure of the gas 
has also been advanced by Mac Low et al. (2005) and 
Avillez & Breitschwerdt (2005). All such simulations 
have ignored the role of thermal conduction. In this pa- 
per we examine the role of thermal conduction on the 
evolution of the 4- volume of the hot gas in the ISM. This 
question has immense bearing on the phase structure of 
the ISM and has assumed immense topical importance 
in view of the observations and simulations cited above. 
SC92 and SC93 used a simplified representation of 
the magnetic field in their early calculations of the evo- 
lution of SNRs with thermal conduction, and wc im- 
prove on that here. Velazquez et al. (2004) include the 
anisotropy in the thermal conduction introduced by the 
magnetic field, but do not include the Lorentz force of 
the magnetic field on the gas. The inclusion of ten- 
sion forces from magnetic fields makes the evolution of 
the SNR anisotropic (Balsara et al. 2001). Furthermore, 
in the presence of a magnetic field, thermal conduction 
itself becomes anisotropic. Thus, the geometry of the 




Fig. 1. — Log density overlaid with velocity vectors (a 
and c) and log temperature overlaid with magnetic field 
lines (b and d) for Run LI at 60 kyrs. The upper pan- 
els exclude thermal conduction, while the lower panels 
include thermal conduction. The tick marks are 10 pc 
apart. 

magnetic field cannot be neglected when calculating the 
thermal heat flux in a multidimensional calculation. 

In this Letter we report on the results of a series 
of MHD simulations on the evolution of SNRs with 
anisotropic thermal conduction. We explore a param- 
eter space of initial ISM conditions in density, tempera- 
ture, and magnetic field strength (Section 2). In Section 
3 we analyze the effects the environment has on the evo- 
lution of the remnant. We examine the porosity of hot 
gas in these different environments in Section 4. Section 
5 provides some conclusions. 

2. Numerical Setup 

Our simulations are performed using the RIEMANN 
code; see Balsara (2004) and references therein. We dis- 
cuss the details of our treatment of anisotropic thermal 
conduction with tests and applications to SNRs in an ac- 
companying paper (Balsara et al., in preparation). We 
draw on Cowie & McKee (1977); Balbus (1986); SC92 for 
our formulation of the classical and saturated thermal 
conduction fluxes. Wc utilize a flux limiter to transition 
from the classical flux to the saturated flux (SC92). The 
transition from classical to saturated conduction causes 
the flux to change character from parabolic to hyper- 
bolic (Balbus 1986). We account for this with a Newton- 
Krylov multigrid method that correctly incorporates the 
anisotropic structure of the thermal conduction operator 



in the presence of a magnetic field. 

We solve the MHD equations in cylindrical geometry, 
using a mesh that is 384x384 uniformly spaced zones in 
the R-z plane. We use the cooling function of MacDon- 
ald & Bailey (1981) for radiative cooling, with heating 
balancing cooling in the quiescent ISM. 

We choose our initial conditions to span a range of 
interesting ISM densities and temperatures. The values 
of the density, temperature and magnetic field in the 
initially quiescent medium are given in Table 1. Runs 
LO and HO label SNRs propagating into low and high 
density ISMs without magnetic fields. Runs LI and HI 
catalogue similar propagation into ISMs with a magnetic 
field of 3 fj,G. We initialize our SNRs with 10^"'^ ergs of 
energy. The physical size of our grid is chosen to capture 
the outer shock through the duration of our simulations. 

3. The Structure of SNRs with Thermal Con- 
duction 

Fig. 1(a) and 1(b) show the logarithmic density with 
velocity vectors overlaid, and the logarithmic tempera- 
ture with magnetic fields overlaid, for Run LI without 
thermal conduction, at a time of 60 kyr. Fig. 1(c) and 
1(d) show the same for Run LI, now repeated with ther- 
mal conduction. We notice that the location of the outer 
shock in the radial direction is unaffected by the inclu- 
sion of thermal conduction. We also observe that at this 
early epoch the hot gas bubble with thermal conduction 
is somewhat smaller in Fig. 1(c) and 1(d), showing that 
thermal conduction does reduce the size of the hot gas 
bubble. This difference becomes even more exaggerated 
at later times. We also observe that the mean density 
in the runs with thermal conduction is higher, and the 
mean temperature is correspondingly lower. Fig. 1(b) 
shows that the hot gas bubble has a temperature up to 
10* K, while Fig. 1(d) shows that the temperature in 
the bubble when thermal conduction is included is at 
most lO'^ K. Thus we see that the inclusion of thermal 
conduction predisposes the hot gas bubble to be strongly 
radiative in soft x-rays, a feature not found in the sim- 
ulations that exclude thermal conduction. 

We show the time evolution of the mean density and 
temperature of the hot bubble in Fig. 2, along with 
a comparison with the evolution of density and tem- 
perature for the same initial conditions, but without 
thermal conduction. The mean densities and temper- 
atures shown in Fig. 2 include all gas with tempera- 
ture T > 10^ K. We see that the inclusion of thermal 
conduction produces a half-order of magnitude greater 
density in each of the SNRs studied here. Likewise, the 
temperature decreases by half an order of magnitude. 
This is directly related to the fact that thermal con- 
duction exchanges energy between the colder and hotter 
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Fig. 2. — Time-evolution of the mean density (a) and mass- weighted temperature (b) of the gas in the hot bubble of 
the SNR. This plot compares the results of runs LI and HI with and without thermal conduction present. 
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Fig. 3. — Time-evolution of the volume of hot gas in Run LI and HI, with and without thermal conduction, (a) 
Volume of gas with temperature greater than 7.9 x 10"'' K. (b) Volume of gas with temperature greater than 2.2 x 10^ 
K. 



parts of the remnant. Because the late-time evolution of 
remnants is isobaric, the decrease in the bubble's mean 
temperature is matched by a corresponding increase in 
its mean density when thermal conduction is included. 
The curves terminate after ~ 2 Myr (in the high-density 
ISM) and ~ 7 Myr (in the low-density ISM), as is appar- 
ent from Fig. 2, due to the collapse of the hot bubble. 

4. Four- Volume of Hot Gas 

The amount of hot gas deposited into the ISM con- 
trols its subsequent thermal evolution. As the rem- 
nant matures, cooling and thermal conduction eventu- 
ally cause the hot bubble to cease expanding and be- 
gin to collapse. This is illustrated in Fig. 3, which 
plots the volume of gas greater than 7.9 x 10^ K and 
2.2 x 10® K as a function of time. The relevance of 
these temperature for high-stage ions is explained in the 



next paragraph. The area under these curves provides 
us a measure of the 4- volume at temperatures exceeding 
the above-mentioned values. Comparing the dashed line 
and solid line in Fig. 3(a) shows that the inclusion of 
thermal conduction dramatically reduces the amount of 
hot plasma in the bubble with temperatures exceeding 
7.9 X 10"'' K after 3 x 10^ years. After 1 Myr we see a 
small emergence of hot gas owing to the fact that the 
interior is reheated by shocks in the collapsing bubble. 
Comparing the dotted and dot-dashed lines in Fig. 3(a) 
confirms that a similar scenario prevails for SNRs going 
off in denser media. Fig. 3(b) is even more illustrative 
as it shows that there is a dramatic difference in the 
4-volume of gas hotter than 2.2 x 10® K when thermal 
conduction is included. In this case, the plasma is hardly 
reheated to these temperatures by the shocks driven into 
the collapsing hot gas bubble. 



Table 2: Porosity of hot gas, defined as 5'_i3/V(i)di where S'-is is the Galactic supernova rate in units of 



10 ^^ pc •^ yr ^ . V(t) is the volume of hot gas that exceeds a specified temperature. 
Gas With conduction Without conduction 
Temperature (K) Run LI Run HI Run LI Run HI 




> 2.2 X 10'' 

> 7.9 X IQS 
> 3 X IQS 


5.11 X 10-'*5_i3 3.63 X 10-'^5_i3 
2.09 X 10'^S-i3 2.12 X 10-35_i3 
9.99 X 10-25_^3 3.05 x lO^^s.^g 


4.26 X 10-^5_i3 1-02 X 10" 
9.00 X 10-25_i3 2.66 x 10" 
1.03 X 10-iS'_i3 2.82 X 10- 


-^5_i3 
-^5_i3 
-^5_i3 



The porosity of hot gas is defined in Table 2 and fol- 
lows the formulation of SC93. For situations where our 
runs have similar parameters to those in SC93, we cross- 
checked our measures of the 4-volume with theirs and 
found good consistency. Table 2 shows the 4-volumes 
for different runs for temperatures greater than 3 x 10^ 
K, 7.9 x 10'"' K, and 2.2 x 10^ K. These temperatures 
correspond to the temperatures at which the emission 
is dominated by O vi (measured by FUSE) and O vii 
and O VIII (measured by Chandra). The amount of gas 
injected into the ISM at hotter temperatures is consid- 
erably reduced compared to the 3 x 10^ K gas. Thus, 
the conclusions of SC93, that the filling factor of hot gas 
from SNRs in the ISM is considerably reduced from the 
standard predictions of the M077 model, are strength- 
ened when anisotropic thermal conduction is added to 
models of SNRs. Table 2 shows that thermal conduc- 
tion reduces the gas that could emit in O vi by a small 
factor but the gas that could emit in O vii and O viii 
by large factors. These results translate to protogalactic 
ISMs where they may even infiuence the physics of early 
galaxy formation. 

It has been speculated (Cho 2003; Avillez & Bre- 
itschwerdt 2005; Balsara & Kim 2005) that high rates 
of turbulent diffusion might make thermal conduction 
irrelevant in global studies of the ISM. Balsara & Kim 
(2005) find that very strongly driven turbulence has a 
coherence time of 0.8 Myr. We anticipate Galactic tur- 
bulence to be much more weakly driven, with a larger 
coherence time. Dispersal of the hot gas bubble would 
only take place after several coherence times. Such times 
would be comparable to the lifetime of the hot gas bub- 
ble catalogued in Table 2. As a result, we do not expect 
interstellar turbulence to completely erase the structure 
of the hot gas bubble before it collapses. We therefore 
expect our conclusions to be robust in the presence of 
interstellar turbulence. 

5. Conclusion 

Anisotropic thermal conduction decreases the aver- 
age temperature of hot gas by half an order of magni- 
tude over a period of a few million years, and increase 
the mean density by a similar amount. It also changes 
the 4-volumes of x-ray emitting gas. The emission in 



high-stage ions is also significantly altered by the inclu- 
sion of thermal conduction. This result is important not 
just for simulations of our present Galaxy, but also in 
cosmological simulations of proto-galaxies. 
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